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I. INTRODUCTION 



In the standard model of cosmology, the early universe is described by a homogeneous and 
isotropic Friedmann-Lemaitre-Robertson- Walker model. Small linear metric perturbations 
are responsible for both, the large-scale structure of the present-day universe and the tiny 
deviations from the anisotropy of the cosmic microwave background. The corresponding 
linear perturbation theory has been developed by Lifshitz half a century ago who found 
the solutions when the energy-momentum tensor is that of a perfect fluid. (For recent 
reviews of further developments since see Ref. and 0.) With more complicated forms of 
matter it is typically necessary to resort to numerical integrations of coupled Einstein- 
Boltzmann equations MM- In the case of collisionless matter, some analytic results have 



been obtained by Zakharov and Vishniak but also there, the perturbation equations 
were eventually solved numerically 0,0]. 

In Ref. , a novel framework for the study of cosmological perturbations has been de- 
veloped which is based on thermal field theory. In this formalism the connection between the 
perturbed metric and the perturbed energy-momentum tensor is provided by the (thermal) 
gravitational polarization tensor. Concentrating on post-Planckian and post-inflationary 
epochs, we assume that T <C mpianck and therefore that the gravitational field can be treated 
as a classical gravitational background field. The momentum scale of cosmological pertur- 
bations is set by the inverse Hubble radius ~ T'^/mpianck? which is thus much smaller 
than the temperature. If the particles are furthermore ultrarelativistic, i.e., their masses 
are negligible when compared to temperature, it is only the high-temperature limit of the 
gravitational polarization tensor which is needed to determine the response of the primordial 
plasma to metric perturbations. The leading high-temperature contributions to the gravi- 
tational polarization tensor, which have been calculated first in Ref. (see also Ref. [0), 
describe collisionless ultrarelativistic matter. Using them to provide the right-hand side of 
the perturbed Einstein equations, one obtains self-consistent and manifestly gauge-invariant 



perturbation equations, for which exact, analytic solutions were found in Ref. [p^,[L4Hl6 



In this case, one can show [16| that the perturbation equations are equivalent to a certain 



gauge-invariant reformulation of the Einstein- Vlasov equations |]T^ . 



In Ref. |]T8[ we have started to extend the thermal-field-theory approach to weakly self- 
interacting thermal matter, for which we have chosen scalar particles with quartic self- 
coupling. With the slight generalization to an 0(A^)-symmetric model, the Lagrangian is 
given by 



C{x) = ^-g{x) {\g^^d,(t)d,(t) - - jip^^M'] . (1.1) 

With ^ = 1/6 this Lagrangian is conformally invariant. The precise value of ^ does not 
enter in our calculations, since curvature corrections to the energy-momentum tensor and its 
perturbations are suppressed by a factor T^/rrip^^^^y- <^ 1. We assume that self-interactions 
are much mor important than those, i.e. A ^ ^^/'^pianck- 

The leading self-interaction effects show up as two-loop corrections to the gravitational 
polarization tensor (or, in particle physics terminology, to the thermal graviton self-energy). 
As it is the case with its one-loop high-temperature limit, it satisfies a conformal Ward 
identity, which makes it possible to use momentum-space techniques to evaluate this nonlocal 
object completely by going first to flat space-time and then transforming to the curved 



2 



background geometry of the cosmological model, which in virtually all cases of interest is 
conformally flat. The corrections to the perturbation equations of the collisionless (one- loop) 
case turn out to be such that it is still possible to solve these equations exactly in terms 
of rapidly converging power series. The resulting changes turned out to be perturbative on 
scales comparable to or larger than the Hubble horizon, whereas the large-time behaviour of 
subhorizon-sized perturbations become increasingly sensitive to formally higher order effects 

In the present paper, we complete the derivation of the effects proportional to the scalar 
self-coupling A and go on to include the next-to-next-to-leading terms, which are of order 
A'^/^. To this order, we have calculated the gravitational polarization tensor in Ref. |p!9 



which required the use of a resummed perturbation theory. Besides the necessity to resum 
the induced thermal masses (~ VXT) acquired by soft excitations in the scalar plasma, 
this also requires a resummation of nonlocal graviton-scalar vertices. These effects can be 
included systematically by our quantum-field-theoretic framework, while it is unclear how 
they could be taken into account in a kinetic-theory approach. 

Again, it turned out that the results for the gravitational polarization-tensor which are 
proportional to A^^^ satisfy the conformal Ward identity which allows us to transform them 
to curved space. Although the results are much more complicated than the previous ones at 
order A\ their structure is again such that the perturbation equations can be solved exactly. 
Since the gravitational polarization tensor also satisfies the Ward identities required by 
diffeomorphism invariance, we can continue to work with manifestly gauge invariant variables 
for metric perturbations. Following largely the notation of Bardeen [^, our gauge invariant 



set-up is laid down in Sect. 2 for the background geometry of the radiation-dominated 
spatially fiat Einstein-de Sitter model. In Sect. 3, the solutions for scalar, vector, and tensor 
perturbations are constructed after allowing also for an arbitrary admixture of a perfect 
radiation fluid, with some of the calculational details relegated to the Appendix. As found 
already in Ref. |]TB[, the corrections to the solutions caused by the scalar self-interactions are 



perturbative except when the horizon has grown to become much larger than the wavelength 
(of a Fourier mode) of the perturbation. In Sect. 4 we show that the behaviour of the 
perturbative series can be greatly improved by rewriting it in terms of a Fade approximant. 



This can be tested in the limit ^ oo of our model (|1 . 1| ) . Encouraged by these findings, 
we consider also the late-time behaviour of our solutions in Sect. 5. Sect. 6 summarizes our 
results. 

II. GAUGE-INVARIANT SETUP 

We will consider a radiation dominated Einstein-de Sitter background 

ds^ = S\t) (-dT^ + dijdx'dx^) (2.1) 

with r being the conformal time which measures the size of the horizon in comoving coor- 
dinates. The overall evolution is determined by the cosmic scale factor S{t) which satisfies 
the Friedmann equation 
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The background energy density is related to the pressure by the equation of state p = 3P. 

In flat space-time the ring-resummed pressure of the ultrarelativistic matter described 
by Q reads ^ 



.vr^T^ 15 A 15 / A 

90 \' 8 7^ ^ Y \ ^ 



P = N^{l-^i- + ^{^\ +OiX')\ . (2.3) 



In the Einstein-de Sitter background the same expression holds true, but now with scale 
dependent temperature T{S) = T{S = 1)S~^. The existence of a thermal equilibrium state 
for relativistic matter is guaranteed by a conformal timelike Killing vector field u^/T. The 
corresponding energy momentum tensor 

= P(4n^n, + 6^) , u, = 55j , (2.4) 

is formally that of a perfect fluid and is traceless. This reflects the conformal symmetry of 
the effective action for ultrarelativistic plasmas. 

We make use of the gauge-invariant metric potentials and matter variables introduced by 



Bardeen [gO| . The linear perturbations may be split into scalar, vector, and tensor parts ^ . 
Since we work in spatially flat space-time all variables can be decomposed into plane waves 
with comoving wave number k. Then the linearized Einstein equations may be written in 
terms of the variable 

X := rfc , (2.5) 

which measures the number of (physical) half wave lengths inside the Hubble radius H~^, 
i.e. 

In what follows it will act as both, a dimensionless time variable and a dimensionless measure 
of the size of perturbations. 

To model a two component universe composed of a relativistic plasma (RP) and a perfect 
fluid (PF), we introduce the mixing-factor 

a = . (2.7) 

Prp + Ppf 

Gauge-invariant matter variables X are then decomposed as 

X = aX^p + (1 - a)XpF , (2.8) 
whereas gauge-invariant metric perturbations Y are given by 

y = Frp + . (2.9) 

This split is possible since we shall consider only small linearized perturbations. 
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1. Scalar perturbations 
The scalar or density perturbations obey the equations 

y$ = em (2.10a) 

x'n = 4°^ , (2.10b) 

with em being the density contrast on hypersurfaces that represent everywhere the local 
restframe of matter, and n!^^ is th( 
related to Bardeen's definition by 



restframe of matter, and n!^^ is the anisotropic pressure. $ and 11 are the metric potentials. 



^ 2 

$A = -n - -$ . 

2 

For our purposes another variable for the density contrast turns out to be useful 

6, = e^-it;f . (2.11) 

The matter velocity fp-* (related to the amplitude of the shear |^0[) is given by the solution 
of 

) + l.f = \{em + r^- - - n. (2.12) 

The quantity rj is the entropy perturbation {rj = corresponds to isentropic perturbations). 
It appears as source term in the trace of the perturbed Einstein equations 

a;2 f^" + 1$' + 1$ + -n' - In) = -7] , (2.13a) 

V X 3 X 3 / 

whereas in the 0-0 components we encounter the density contrast eg 

(x^ + 3)$ + 3x$' + en = 3eg . (2.13b) 
Here and in the following a prime denotes the derivative with respect to x. 



2. Vector perturbations 
For vector (rotational) perturbations 




(2.14) 



relates the metric (frame dragging) potential ^ to the matter velocity Vc (relative to the 
normal of constant-r hypersurfaces), which is proportional to the amplitude of the vorticity 
pO| . The anisotropic pressure is given through 

v'. = -l4^- (2.15) 
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3. Tensor perturbations 



Tensor perturbations describe the propagation of gravitational waves, their evolution 
equation has the anisotropic pressure as source term 

{h" + '^H' + H^= 4'^ . (2.16) 



III. PERTURBATIVE THERMAL FIELD THEORY 



A. Fluctuations from thermal A</>^ theory 



To obtain the expressions for the matter variables we use the thermal-field-theoretic ap- 
proach |Tl]] instead of the usual kinetic approach. The perturbation of the energy momentum 
tensor 



7 r^Sgaf^ix ) 

X' Ogaf3{x') 



-9[X) 



n^/^(x,x')5^„/3(x') 



{x)5gaf3ix) 



(3.1) 



is related to the gravitational polarization tensor 



6^r 



x,x ) = 



Sg,^uix)6ga/3ix' 



l5[^-g{x)T^^''{x] 
2 5gai3{x') 



(3.2) 



I.e. 



In the high-temperature domain the effective action is conformally invariant p!9|JT2 
T[S'^g^y\ = r[gfj_iy], therefore the polarization tensor in the conformally flat Einstein-de Sitter 
background reads 



U'"'P''{x,x') 



S 



-2/ 



(2 



s~ 



(3.3) 



Additionally we have performed a Fourier transformation to momentum space, where stan- 
dard thermal-field-theoretic methods apply. Due to the conformal invariance and the in- 
variance under general coordinate transformations, the gravitational polarization tensor has 
only three independent components. We choose them to read 



AiQ) ^ fi°°°°(g)/p, BiQ) ^ n%°(g)/p, ciQ) ^ w\AQ)/p- 



(3.4) 



The explicit expressions for A, B, and C have been calculated in |[T^ to 0{\2). The contri- 
butions through order A consist of "hard thermal loops", i.e., diagrams that are dominated 
by hard loop momenta. These correspond to thermal fluctuations and could have been 
calculated by kinetic theory as well (despite the collision term that cannot be derived by 

3 

flrst principles). The order A 2 stems from the necessary resummation a la Braaten-Pisarski 
||2^ where propagators and vertices of the scalar particles have to be dressed in order to 
avoid infra-red singularities. The resummation involves a particular (inflnite) subclass of 
Feynman diagrams of ordinary perturbation theory which presumably lie far beyond the 
scope of kinetic theory. A, B, C and their Fourier transforms are listed in Appendix A. 
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1. Scalar perturbations 



From and the definition of tlie gauge-invariant matter variables pO| it follows that 



^RP = 



1, 



e, Rp = -2$ + 4^[A-^]*($ + n), 



,(0) 



i,p^3z^[^(A--)]*($ + n) 



(0) 

T RP 



-i8^[(a.2--)(A--) + -]*(<i. + n) 



(3.5a) 
(3.5b) 

(3.5c) 

(3.5d) 



The * denotes the convolution {g * f){x) = dx'g{x — x')f{x') and the operator defines 
the Fourier transformation 



J-'[g]{x) = lim — / due' 

7^0+ ZTT J-oo+i-y 



(3.6) 



The particular choice of the integration contour corresponds to retarded boundary condi- 
tions. 

After performing the Fourier transformations (Appendix A) the matter variables are 
written in terms of the integral kernel 



7^(0) ( 



X] 



MX) + 



5A 



{2k - jo - cos) (x) 



+ ir (4 I ( cos+-J_i + 2 Jo + -Ji + - -3o 



77^ 



X 



-2v" - 4 



U — K 



K 



(3.7) 



according to 



eg RP = 2$ + 4n - 4is:(°) * ($ + n)' 

,(0) -37^(0)'* ($ + n)' 



TT. 



-'s RP 
(0) 



,^ ^ -6(ir(°) + 3ir(°)") * ($ + n)' 

At the origin this kernel behaves like 

ir(°)(o+) = i, kW'(o+) = o. 



(3.8a) 
(3.8b) 
(3.8c) 

(3.9) 

The scalar matter variables (|3.5| ) couple only via the sum of the metric potentials $ + 
n = $Ar to the nontrivial component A of the polarization tensor. $Ar is the only scalar 
contribution to the electric part of the Weyl curvature tensor |^3[ and can be interpreted as 
the generalization of the Newton gravity potential. There is no scalar contribution in the 
magnetic part of the Weyl tensor. Were it not for the conformal symmetry of the effective 
action, the additional components of the polarization tensor would also couple to 11, the 
potential for the anisotropic pressure. 
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For vanishing entropy perturbations, Eq. ( p.l3a| ) can be integrated for $, 
2 r^' 

= / dx' (x' cos(x - x') + 4sin(x - x')) $'/v(x') 

3x J 



^ ix cos(x — x) + sin(x — x')) $Ar(x') 

3x 



+ + (3.10) 



X 



X 



where Ci,C2 are integration constants. Together with Eqs. ( p. 10b ) and ( 3.8c ) we can also 
obtain a single equation determining $Ar (see Eq. (|3.21| )). 



2. Vector perturbations 



In a similar way 



Vc RP = + \:F[{uJ^ - 1)(A - i) - 5] * VP 



(3.11) 



is derived. With 

ir«(x) = 



1,^ 15A 15A3/2, 

77 1 + 



15 / A \ ^ /I 



(3.12) 



and 

irW(0+) = --, K«'(0+) = 
3 

the velocity amplitude is given by 

7;eRP = -3(ir«')*^- 



(3.13) 



(3.14) 



3. Tensor perturbations 
The evaluation of ( p.l| ) for the anisotropic pressure leads to 

4\v = - -\)- 4(c.2 - 1)5 + 2C - + 3] * 



(3.15) 



with the kernel 

K^^\x) = 



j2{x) ^ 15A _ 15A3^^ ^ 5A ji(a;) 



x^ 



27r3 



X 



+ J0+J2][X 



. (3.16) 
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B. Initial conditions 



To specify the dynamical equations for cosmological perturbations completely, initial 
conditions for the metric potentials are needed. We fix them in the limit Xq 0, because 
this will allow us to derive exact solutions of the integro-differential equations in terms of 
generalized power series [0. The initial conditions are given by a set of numbers 7^'^-'. 
These are related to the n-th momenta of the particle's distribution function |[l^. Here they 
show up in the convolution integrals because the lower integration boundary xq truncates 
the support of the metric potentials to a half space. This is done in the convolution integrals 
by the replacement 



00 

Y'{x) Y'{x)e{x -Xo) + J2 In^^^^'H 

X — Xq ) 



n=0 



(3.17) 



where Y stands for $7v, \f , and H, respectively. 

As is well known from simple forms of matter like perfect fluids [|1| and from work by 
Zakharov and Vishniac on collisionless matter, there are two branches of solutions, 
regular and singular ones. This is also the case for self-interacting plasmas. We will con- 
centrate on regular solutions in our evaluations, but will nevertheless sketch the behaviour 
of the singular solutions below. The singular solutions necessarily violate the assumption of 
a Friedmannian singularity, but they are relevant when fitting the evolution of cosmological 
perturbations to a previous epoch at some small nonvanishing value of r p3,|2H]. A detailed 



discussion of singular solutions for collisionless matter is given in |16]. Due to geometrical 



effects |ro|], the singular solutions permit superhorizon oscillations. This effect depends on 
the ratio a only, and is not at all sensitive to the 7^°-' . Only the normalization of the singular 
part of the solutions has to be fixed by an initial condition. 

The regular solutions are determined by the 7^"^^ We will restrict our attention to 
isentropic (adiabatic) perturbations, which may be left over from an earlier inflationary 
epoch. We will neglect all 7^**^ with n > 2. This is motivated by the findings of |TB| that 
these quantities are related to the higher moments of the kinetic distribution function. The 
restriction to n < 2 means that we fix the momenta which directly occur in the energy- 
momentum tensor and set all others to zero initially. 



1. Scalar perturbations 
From ( ^.101 ) and ( |2.12| ) the small x behaviour of the matter variables 



,(0) 



x^^O) 

x($(o) + 2n(o)) 

x^U{0) 



follows. Therefore, vf'\0~^) and 'it!j^\o^) = for regular solutions. This means that the 
r.h.s. of ( |3.8b|) and ( p.8c|) after performing (|3.17| ) have to vanish as well, i.e.. 
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g^{0);^(0)(n+l)^Q+) (3_^g) 
n=0 

oo 

J2 iTiK^""^ + 3ir(°)")^"H0+) = . (3.19) 

?i=0 

Due to 

ir(°)'(0+) = + 3ir(°)")(0+) = (3.20) 

we are free to choose any 70°''. It is remarkable that these relations hold for collisionless 
matter and are not changed by the addition of weak self-interactions. A completely arbitrary 
choice of the 7^°-' is not possible. In the following we choose 70°^ 7^ and all higher 7^°) 
vanishing. This corresponds to the usual choice for isentropic initial conditions, e.g. [EB . 



In the above mentioned equation for <l>7v the arbitrary constants Ci and C2 remain to be 
fixed. Initial conditions are included by the replacement ( ^.17[ ) which yields, together with 
Eqs. ( p.lO| , |2.10E| , p.8c| ), the integral equation for 

$Ar(x) = -- / fix' (x'cos(x - x') + 4sin(x - x')) <l>'^(x') 

r c/x'(K(°) + 3K(°)")(x - x)^'i,{x') + i.e. (3.21) 
Jo 



X Jo 
18a 



x"^ Jo 



with the initial conditions 



i.e. = 3Ci - 2<|.^(0) + 4 Y: 72n(-)"(n - 2) + ^ ( 3^2 + 2 ^ 72n+i(-)"(2n - 3) 



X 



n=0 / \ n=0 



1 Rn/ 

J2 ln{K^'^ + 3ir(o)")(x))("). (3.22) 



n=0 



begin For regular solutions, C2 is determined by requiring that the term in the second round 
brackets vanishes, and Ci is related to the initial value $Ar(0) by the constant term in the 
series in x. 

The metric perturbation $rp is given through 

$RP = 5 {e, KP + ^v^X) ■ (3.23) 

The regularity requirement establishes a relation between $jv(0),n(0) and the 72°'' coeffi- 
cients. Using this relation together with the expressions for $Ar(0) and 11(0), Ci may be 
expressed as a sum of coefficients 72°'' (see Appendix B). 



2. Vector perturbations 

From ( 2.15|) a constant Vc pf follows, which is a consequence of the Kelvin- Helmholz 
theorem. For a perfect fluid alone, this forbids regular solutions at all, because it entails 
a singular frame dragging potential This situation changes when collisionless or weakly 
self-interacting matter is added (a discussion of potential cosmological consequences like 
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primordial magnetic fields can be found in ^7j). One can have a regular vector perturbation 
sustained by the relativistic plasma alone, or one can compensate a primordial vorticity 
in the perfect fluid component, which has constant Vc pf, by a nonvanishing vorticity of 
opposite sign in the relativistic plasma, and a growing net vorticity is generated by the 
nontrivial evolution of the latter. 
Eq. (|2.14| ) implies for small x 

Vc = avc Rp + (1 - a)vc PF ~ a;^ . (3.24) 

Therefore, ( glil) leads with ( CT ) to 



..PP = -^E7i^^i^^^^^"H0^). (3.25) 

1 - « n=0 

As can be seen from the above formula, in the absence of a primordial perfect fluid vorticity 
one has to have nonvanishing coefficients 7^^-*, n > 1 for nontrivial regular solutions. In 
the presence of a perfect fluid component, we shall restrict ourselves to nonzero 70^^ and 
vanishing higher coefficients. 

3. Tensor perturbations 
Eqs. (|2l6|) and ( glSD with (|3T^ ) lead to 



^7(2) K{2)W(o+) = . (3.26) 



n=0 



Because tensor perturbations correspond to gravitational waves, there are nontrivial solu- 
tions for all 7^^^ vanishing. Then the initial conditions are specified by the amplitude H(0) 
and its derivative H'{0). This will be our choice in what follows. 

C. Solutions through 0{X^/'^) 

The singular solutions to the equations for the dynamics of cosmological perturbations 
are obtained with the Ansatz 

yW(a;) = x'"F('*)(x) (3.27) 
with F(")(0+) ^ and finite (a = 0, 1, 2). This yields 



5 1 I a 



Therefore, for a greater than 
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a takes complex values, and thus gives rise to superhorizon oscillations ~ cos([Im(T] Inx). 
For small A > the value of ctcrit. is increased. This is as one may expect, since in the 
collision-dominated case of a perfect fluid superhorizon oscillations do not occur. For large 
values of A, A > 87r^/63 ~ 1.25 this trend is reversed, but there the perturbative series can 
no longer be trusted, because the A^^^ correction begins to dominate over the A^-term. 

Turning now to the regular solutions, the power series Ansatz detailled in Appendix B is 
made. This leads to recursion relations for its coefficients that can be solved as also given in 
the Appendix. These have been evaluated with a Mathematica code On superhorizon 
scales the effects of the self-interactions on the regular solutions are small. E.g., for scalar 
perturbations with a = 1 and the above initial conditions, 

-?\0)^-|fl-iA^I^fA| ^0(,.,| (3.30) 



e 



(0) ' ^ 7 28 7r2 16 I tt^ 



m 



the amplitude of the anisotropic pressure becomes smaller when the plasma is collisional 
for values of A smaller (47r/21)^ ~ 0.36. Since for tightly coupled plasmas we expect the 
behaviour of a perfect fluid, i.e. the above ratio should be zero, we conclude that the 
perturbative result through order A'^''^ can be trusted only for A < 0.36. 

Some solutions for the density contrast in the subhorizon region are shown in Fig. |l|. 
The three lines show the behaviour of \em Rp(a;)| for an ultrarelativistic plasma alone. We 
plotted solutions for A = 0,1/2, and 1. The coUisionless solution (full line) decays due to 



directional dispersion p9[. The phase velocity of the oscillations is the speed of light. For 
very small values of A no visible effect occurs in the plotted region, whereas for bigger values 
a considerable change in the damping behaviour and in the phase velocity is observed. For 
A = 1 (dotted line) the density contrast starts to grow again at a;/7r ~ 6, which in fact is 
associated with rather unlikely beats. Also for the smaller value of A = 1/2 such a behaviour 
arises for large enough x. This will be analysed later in the next section. 
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em Rp(a;) 



. 5 

0.2 
. 1 
. 05 

. 02 
.01 




FIG. 1. The full line shows the density contrast \€m rp| for collisionless matter in the subhori- 
zon region. The dashed and dotted solutions show the effect of coUisional matter through 0(A^/^) 
with A = 1/2 and 1 respectively. 

IV. FADE IMPROVEMENT 

As we have seen, in order to have sizable effects from the self-interactions of the thermal 
matter, we have to adopt sizable values of A. However, the perturbative results quickly 
become unreliable with increasing A. In order to get a somewhat more quantitative idea of 
the problem and of ways to improve this situation, we first inspect a solvable case. 



A. The thermal mass in the limit N 



oo 



It is well-known ||3^ that in the limit ^ oo the model of Eq. ( |1 . 1|) becomes exactly 
solvable. In this limit, the thermal mass of the scalars is exactly given by the resummed 
one-loop gap equation 



6A 



dx 



TC^ Jm/T — 1 

The first two terms of the perturbative series 



m' 



•/T^ = X--X^/^ + 0{X^\nX) 

71 



(4.1) 



(4.2) 
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are naturally a good approximation for very small A, but apparently breaks down when 
A ~ 1. Indeed, for A = 1, Eq. (|4.2| ) gives m? /T"^ k, 0.05, which is an order of magnitude too 
small when compared with the result following from Eq. ( ^.11 ), which yields m^/T^ 0.53. 
However, rewriting Eq. (|4.2| ) in the perturbatively equivalent way of 



A 



1 + ^AV2 



+ 0(A2lnA) 



(4.3) 



considerably extends the range of A over which the first two terms of the perturbative series 
give a faithful picture of the actual behaviour. For A = 1, Eq. ([4.3|) yields m^/T^ ^ 0.51, 
which is only a few percent too small. 




A 



FIG. 2. m? /T"^ as a function of A in the exactly solvable case of — > oo (full line). The 
long-dashed curve is the perturbative result (4.2); the short-dashed one the Fade-improved (4.3). 



Of course, for still larger values of A, also Eq. ( [4.3|) becomes increasingly imprecise, but 
the improvement over Eq. ([4.2| ) remains striking. What we have done by going from Eq. ( |4.2| ) 
to Eq. ([4.3| ) can be viewed as replacing the first terms of a power series in a/A 

X = a + 6A^/2 + cA + rfA^/2 + . . . (4.4) 

by its (2,1)-Pade approximant (again in powers of \/~\) ||3T| 

In the case of the thermal mass, a = since we have no tree-level mass to start with, and 
6 = since the plasmon effect is down by one-half order with respect to the perturbative 
one-loop one. Extending this procedure to other quantities, we thus still have 6 = 0, so that 
the results through order A^^^ determine the four parameters of the corresponding (2,1)-Pade 
approximant s, 

a = a, P = —ad/c, 7 = c, 6 = —d/c. (4.6) 
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B. The resummed 1-loop kernel 



In our applications, the accuracy of the perturbative result will not only deteriorate when 
A is increased, but also when x becomes too large, as we have seen in Fig. [1|. This is so 
because the asymptotic behaviour of the Fourier transforms (|3.7|J3.12| , p.l6| ) with a; — > oo is 



such that for any small but finite value of A the correction of order A^^^ in K'^"''> eventually 
overtakes the A^ term which in turn overtakes the lowest order term. Apparently, this signals 
a breakdown of perturbation theory at large x. 

The origin of the difference in the asymptotic behaviour in x comes from an increasingly 
singular behaviour of the discontinuity of the functions A, B, and C at higher orders in A. 
For example, from ([Al|) one notices that at lowest order the discontinuity of A across the 
branch cut between = ±1 is a constant; at order A^ it is logarithmically singular at = ±1 
and in addition, there are now simple poles at these points; at order A'^''^ the singularities 
are still worse. Since these singularities occur at the end points of the integration region 
contributing to the Fourier transform, they become dominant for the large x behaviour of 
the latter. 

On the other hand, these singularities at the light-cone should not exist at all since 
the originally massless scalars have acquired thermal masses. Indeed, keeping the thermal 
masses in the integrals without expanding them on account of them being proportional to 
A, shows that the complete discontinuity is smooth at u = ±1 and that also the simple 
poles there are spurious [|l8l. For instance, the one-loop contribution to K^^\x) = jo{x) is 



modified by including the thermal mass in the scalar propagators according to 

(4.7) 




up to terms whose amplitude is suppressed by explicit powers of A. With a non-zero m = 
\/AT, the integrand is now seen to vanish at oo = ±1. Instead of being a constant, it vanishes 
at the endpoints of the branch cut together with all its derivatives, but rapidly recovers the 
bare one-loop value away from the light-cone. This is a negligible effect for small x, but the 
large x behavior is changed completely. 

Eq. ( |4.7| ) can be evaluated by a Mellin transform which yields 

Krix) = X (4.8) 

V a/2 (-l)^+^(fc-4)C(fc-3) /2yi-'=)/2 

to {2nr-ni + k/2) [xJ ^(i-'^V^l^)' 

where in the term with k = 4 one has to substitute {k — 4)C(/c — 3) — >• 1. For small x, 

Krix) = (^Ux) - ^ cos(x) j + 0(A^/2) (4.9) 

is a good approximation; for x ^ 1, on the other hand, the complete function Kl'^^'{x) turns 
out to decay even faster than jo{x), oscillating with a reduced phase velocity 
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(4.10) 



A better approximation is thus obtained by modifying 

This indeed resums and thus removes the terms of order A in K^^'^ which are the most domi- 
nant for large x, leaving only such which are logarithmically larger than Ki'^ However, the 
kernel K^^^ has to satisfy (|3.20|) in order that the equation for scalar cosmological perturba- 



tions remains integrable. In Ref. |]18[ we have found that this can be fulfilled by adjusting 
the order A term in K^^^ so that it contains the same phase velocity as the modified k[^^ 
and correcting the prefactors without modifying the lowest orders in A. 

However, it is difficult to see how the effect of the higher terms in Eq. (|4.8| ) can be 
accounted for by similarly simple modifications. The functions that come with higher powers 
of A are in fact more and more diverging with x and only the infinite sum is a decaying and 
regularly oscillating function — any truncation results in a function that oscillates with 
interferences before exploding eventually. 

Unfortunately, we are able to compute only the first few terms of the contributions to the 
kernel, Eq. ( [4.7| ) being only one particularly simple contribution. In order to improve the 
behaviour of our results at large values of x, we propose to use the same Fade approximation 
that was working rather well above. The most obvious way for doing that would be to 
turn the coefficients in Eq. (|4.6| ) into functions of x. However, this turns out to violate 
the integrability constraints on K^"'^ A method which manifestly respects the latter is to 
perform a Fade improvement on each of the Taylor coefficients of K^""^ as a function of x. 

In order to test this procedure, we have applied it to the first three terms of the expanded 
kernel ( |4.8| ) . With A = 1 , the result is shown in Fig. ^. There the full line gives the complete 
function (|4.7] ) and the long-dashed line shows the truncated series through A^^^. The latter 
deviates quickly from the complete function and after few oscillations is completely off. The 
Fade improved function, where each term in the power series in x is replaced by a (2,1)- 
Fade approximant, behaves instead much more regularly and is quite close to the complete 
function for all values of x. 
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FIG. 3. The function (full line) for A = 1 and two perturbative approximations: The 

long-dashed line gives the first three terms of its perturbative expansion; the short-dashed one its 
Fade improved version. 

We therefore expect that the analogous procedure for improving the kernels we have 
obtained perturbatively up to and including order A^''^ allows us to extend the range in both, 
A and x, where the solutions to the corresponding equations for cosmological perturbations 
can be trusted, to about A ^ 1 and x/tt ^ 10. 

V. PADE-IMPROVED SOLUTIONS 



On superhorizon scales the perturbative results became obviously unreliable already at 
moderate values of A because, among others, the 0(A^/^)-contributions in Eq. ( |3.29D stopped 
acrit. from increasing beyond A ~ 1.25. With the Pade-improved results we have 



acrit.(A) 



4 + 21^ 



2^4 + 21^-54 



(5.1) 



which shows an ever-increasing behaviour up to very large A. Likewise the Pade-improved 
version of Eq. ( |3.3UD is a monotoneous function. In both cases, the results for the weakly- 
interacting plasma remains far from the perfect-fluid ones for A ~ 1; only for A ^ 10^ these 
would be reached, where a perturbative treatment is certainly inadequate. 

In the following we shall inspect our Pade-improved solutions for A = 1 and x/tt > 10. 



17 



A. Scalar perturbations 



In Fig. ^, the density contrast associated with a scalar perturbation in a pure relativistic 
plasma is given for the interacting and the coUisionless case. The difference turns out to 
be moderate so that a perturbative treatment seems justified. The main effect turns out 
to be a somewhat decreased phase velocity and a somewhat diminished exponent in the 
power-law decay. This is exactly what one would expect in view of the behaviour of density 
perturbations in the presence of a perfect fluid, where the phase velocity equals l/\/3 and 
where damping through directional dispersion is inoperative. 




FIG. 4. The density contrast is shown for coUisionless matter (full line) and for an ultrarela- 
tivistic plasma with A = 1 (dotted line), which is the Fade-improved solution. 

In Fig. ^ the density perturbations are shown for a two-component system with an equal 
amount of perfect fluid and relativistic plasma. There is little difference from the case 
considered in Ref. |T^, where the relativistic plasma was coUisionless. The main effect is 



again a diminished phase velocity, which is exhibited in the magnified Fig. y. There is no 
longer a simple decay-law for the density perturbations in the plasma component, because 
it is strongly influenced by the comparatively large over- and underdensities created by the 
acoustic waves propagating in the perfect fluid component. 
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FIG. 5. Density perturbations for a mixture of a perfect fluid (dashed line) and an ultrarela- 
tivistic plasma (dotted line) with a = 1/2 and A = 1. 
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FIG. 6. For the same mixture as in Fig. 5 the subhorizon perturbations of the plasma com- 
ponent are shown for A = (full line) and A = 1 (dashed line). 
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In Fig. 0, the anisotropic pressure associated with the scalar perturbations in the two- 
component case is given, compared with the colhsionless version. Except for the third peak, 
there is remarkably little difference between A = and A = 1, although one might have ex- 
pected that the anisotropic pressure would be the most sensitive quantity to self-interactions 
in the plasma, since a collision-dominated perfect fluid forbids anisotropic pressure com- 
pletely. 



2 




FIG. 7. The anisotropic pressure is plotted for a = 1/2. The full line shows the behaviour for 
a plasma with A = 0, whereas the dotted line shows the solution for A = 1. 



B. Vector perturbations 



As we have mentioned, the presence of a relativistic plasma opens the possibility of having 
regular vector perturbations. In a perfect-fluid alone, the velocity amplitude is constant in 
the radiation-dominated epoch by virtue of the Kelvin- Helmholtz circulation theorem 



which leads to a singular behaviour of the frame dragging potential. Adding in a plasma 
component with initially compensating vorticity, however, allows nontrivial regular solutions. 
These solutions were first found within the thermal-field-theoretical treatment and it was 
pointed out in Ref. [27| that they might have interesting applications in the open issue of 
primordial magnetic fields. 

In Fig. ^ such a solution exhibiting the generation of a net vorticity which approaches 
a constant velocity amplitude is given for A = 1. In this case there is even only a relatively 
small deviation from the colhsionless scenario without the Fade improvement (Fig. |^). With 
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it, however, the difference becomes even rather tiny. This is also somewhat unexpected, 
since the very existence of these sohitions hinges on having a plasma component that is 
approximately collisionless. The effect of self-interactions are found to give only a small 
increase of the period of the wiggles in the vorticity of the plasma component while it dies 
from directional dispersions. 



1 



0.5- 




0.1 0.2 0.5 1 2 5 10. 



FIG. 8. In a universe containing a mixture (a = 1/2) of a perfect fluid and an ultrarelativistic 
plasma rotational perturbations li^c(a;)| may survive in the subhorizon region (full line). The dotted 
line shows the rotational perturbation of the plasma component \vc kp{x)\ with A = 1. 
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FIG. 9. For the same mixture as in Fig. 8, the rotational perturbation \vc Rp(a;)| is plotted 
for A = (full line) and for A = 1. In the latter case the dashed line shows the solution through 
order A^/^; the Fade improved solution is shown by the dotted line. 

Without a perfect-fluid component, one can have a regular solution which has a growing 
velocity amplitude on superhorizon scales and which decays from directional dispersion after 



horizon crossing, which is shown in Fig. |10| In Fig. |TT], a magnified picture of the subhorizon 
behaviour is given, which shows both a small decrease of the phase velocity of the oscillations 
and a small reduction of the damping. While this is similar to the subhorizon behaviour 
encountered in the scalar case, it could hardly be anticipated by comparison with the perfect 
fluid case, because in the limit of tight couphng these perturbations are forbidden entirely. 
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FIG. 10. Without perfect fluid all rotational perturbations decay on subhorizon scales. The 
full line shows the behaviour of collisionless plasmas (A = 0), whereas the dotted line shows the 
Fade improved solution for A = 1. 
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FIG. 11. The same as in Fig. 10, but just the subhorizon region. 
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C. Tensor perturbations 



Tensor perturbations correspond to primordial gravitational waves. Their large-time 
behaviour is expected to be rather independent of the medium, since it is dictated by energy 
conservation. Indeed, there is only some difference in the behaviour at the time of horizon 
crossing which implies that a relativistic plasma requires stronger initial tensor perturbations 
in order to have equal amplitude in the gravitational waves at late times. There is, however, 
extremely little difference in the behaviour of the solutions for the plasma case for A = 
and A = 1, see Fig. The self-interacting case is closer to the perfect-fluid case, but only 
very little so. 

As Fig. |l^ shows, the perturbative result is moreover rather insensitive to the Pade- 
improvement. 
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FIG. 12. The amplitude of a gravitational wave for a a perfect fluid (dashed line) and 

an ultrarelativistic plasma. The full line solution corresponds to A = and the dotted line to the 
Fade improved solution for A = 1. 
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x/tt 

FIG. 13. As in Fig. 12, but only for a = 1 with A = 1. The full line shows the solution through 
order A^/^ and the dotted line the Fade improved solution. 

VI. CONCLUSION 

We have studied the effects of weak self- interact ions in an ultrarelativistic plasma on 
cosmological perturbations through order A^^^ in A0^-theory. At this order it turns out 
that perturbation theory requires a resummation of an infinite set of higher-order diagrams. 
This is natural to include in the thermal-field-theory approach to cosmological perturbations, 
which in the coUisionless case is equivalent to the usual approach based on classical kinetic 
theory, but now leaves the latter clearly behind. 

While it still turned out to be possible to exactly solve the perturbation equations by 
means of a power series ansatz, we have found that the relatively large coefficients of the 
order- A^/^ corrections make the results reliable only for rather small values of A, and even 
then there is a breakdown of perturbation theory in the asymptotic late-time behaviour. The 
latter comes from increasingly singular contributions to the gravitational polarization tensor 
for light-like momenta and could be cured by a further resummation procedure. However, 
the effects of the latter turned out to be well approximated by a (2,l)-Pade-improvement 
of the perturbative results, which also drastically improves the apparent convergence of the 
results for smaller times. 

The concrete results obtained showed a tendency toward perfect-fluid behaviour, but 
with A = 1 all of them are still (perhaps surprisingly) close to the coUisionless case. The 
main effects turned out to be a small increase of the critical mixing factor of a perfect-fluid 
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component with a relativistic plasma above which one can have singular solutions exhibiting 
superhorizon oscillations. Concerning the regular solutions, we have found a decrease of the 
phase velocity of scalar perturbations and a reduction of its damping. In the case of vector 
perturbations (corresponding to large-scale vorticity), which are only possible in the presence 
of a relativistic plasma, similar effects were found, but quantitatively much smaller. Very 
little effects from self-interactions were finally observed in the case of tensor perturbations, 
which correspond to primordial gravitational waves. 

From this one may conclude that a description of a primordial plasma built from weakly 
interacting elementary particles through perfect-fluid models is in general applicable only for 
scales far below the Hubble radius. At the scale of the horizon and beyond, self-interactions 
can be treated perturbatively, at least in the model considered here, and there can be 
significant differences from the perfect-fiuid behaviour, in particular in the case of rotational 
perturbations. 
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APPENDIX A: GRAVITATIONAL POLARISATION TENSOR COMPONENTS 

AND ITS FOURIER TRANSFORMS 



The Fourier transform of the combinations of the gravitational polarisation tensor com- 
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with powers of a;, as required in 3.11 , 3. 15 ), boils down to the evaluation of the following 
integrals 
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where additional powers of uj in the integrands can be obtained by suitably differentiating 
those functions with respect to x. Let us first consider a typical integral appearing in the 
power-series expansion in x of v{x), 
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By integration by parts with respect to the root, it can be shown to satisfy the recursion 
Sm = 2m{Sm-i — Sm) + '^V^^^i^ + TTi) / ml which may be solved for Sm, 

^^-2V7r2^^,^^^2m-2/)' 



that finally leads to the expression 
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can be removed by integrating by parts in the same manner as above with the remaining 
integral being of the type of Sm with the root in the numerator. satisfies the recursion 
3Tm = (2m + l)(Tm_i — Tm) + 4(5*^ — Sm+2) which can be solved using the initial integral 
To = — 2n/3. Splitting off the contributions oc tt^ which can be reexpressed in terms 
of Ji{x) finally yields 

^(x) - 6 V (-ir^^'" + / - r(/ + i) J. r(j + i) \ 

^^^^ (2m)! r(i)r(m + 3)l' fro 4r(/ + f) j!(/ - J + |) j ' ^ ^ 



APPENDIX B: RECURSION RELATIONS 



In order to solve the integro-differential equations for the metric potentials a power series 
ansatz is appropriate. The integral kernels are expanded as 
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In this Appendix we restrict our attention to regular even solutions, as the kernels are even 
functions of x. This means that we assume 72^+1 = for all n and a. Regular odd and 



singular solutions may be obtained in a similar way ||16|| . 
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1. Scalar Perturbations 



With 
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X 



2n 



n=0 



(2n)! 



Eq. ( pj.21D leads to the recursion relation 
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for n > 1. The regularity at x = demands 
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Accordingly, the metric perturbation 
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is determined from Eqs. (|2.10b| ) and (|S.8c|) . Its recursion reads: 
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for > 0. The relativistic plasma component perturbation 
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follows from ( ^■2S ) and reads 
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for n > 0. From regularity at x = 
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follows. 
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2. Vector Perturbations 



The vector metric potential 



X 



2n 



is determined by the Eqs. ( p.l4|) and (|3.14| ). The corresponding recursion is 
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for n > 0. The relation ( |3.25|) couples the vorticity of the perfect fluid component to the 
coefficients 72}''. 



3. Tensor Perturbations 



For the metric potential H the ansatz 
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is made. Regular even solutions follow from Eqs. ( p.l6|) and (|3.15| ) by 
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for n>l. The initial value ho specifies together with the inhomogeneous terms the regular 
solution. 
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